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Abstract 

Within the chiral constituent quark model of Glozman and Riska, we discuss 
the stability of heavy pentaquarks, i.e. hadrons containing four light quarks 
and a heavy antiquark. The spin-dependent part of the Hamiltonian is dom- 
inated by the short-range part of the Goldstone-boson-exchange interaction. 
We find that these systems are not bound, having an energy above the lowest 
dissociation threshold into a baryon and a meson. 

The question of the existence of hadrons other than mesons qq and baryons g 3 remains 
open. Multiquark states as tetraquarks q 2 q 2 , pentaquarks q 4 q or hexaquarks g 6 , could exist in 



* Supported by the EU program ERBFMBICT 950427 



1 



QCD. Based on a model with a chromomagnetic interaction, Jaffe |J predicted the existence 
of several multiquark states and among them the uuddss hexaquark, with J p = + and 
1 = 0, known as the if -particle. Since then, production and decay properties of such exotic 
hadrons have been derived within several models and extensive efforts have been made to 
find them experimentally 0. The recent high-sensitivity search at Brookhaven || gives no 
evidence for the production of deeply bound if -particles, as predicted by Jaffe. 

This is the present situation in the light sector. However, from theoretical general argu- 
ments [§]-|6|], one expects an increase in the stability of the multiquark system containing a 
heavy-flavour quark Q. Experiments are being planned at Fermilab and CERN to search for 
new heavy hadrons and in particular for doubly charmed tetraquarks [0-0] . Very recently, a 
search for the pentaquark = csuud and P g ~ = csddu production, performed at Fermilab, 
has been reported |Hj . Decay properties into the channel (pirp have been analysed assuming 



a life time ranging from 0.1 to 1 ps and a mass between 2.75 and 2.91 GeV. 

This is an incentive to study the stability under strong decays of pentaquarks q A Q where 
q = u, d or s and Q = c or b. From a theoretical point of view, the field is about ten years old. 
In the constituent quark model based on one-gluon exchange (OGE) interaction between 
quarks, the existence of a stable compact pentaquark was first predicted in Refs. ]TT|jr2 



provided that at least one of the light quarks is strange and Q = c or b. But with a proper 
treatment of the kinetic energy and a realistic breaking of the flavour SUp{3) symmetry, 



the pentaquark turned out to be unstable [13,14]. A systematic study of all pentaquarks 



has been performed in Ref. [15] where over a dozen are found to be candidates for stability. 



Among them those with strangeness S = — 1 or —2 were the most favourable. Calculations 



using an instanton [16] or a Skyrme model [17] indicate that pentaquarks, which are not 



necessarily strange, can appear as bound states or as near-threshold resonances, depending 
on the parameters of the model. On the other hand, the possible existence of heavy-flavour 
pentaquarks as molecular-type resonances, where the quarks interact by the long-range one- 
pion exchange interaction, has also been considered (Tj|. In those calculations, only the NB* 
system was found to be bound. Finally, a study of charmed multiquarks in a bag model, 



where no bound state is found, has also recently appeared |L9 . 



The present study can be viewed as a natural extension of our previous investigations 
f20|j on the stability of heavy-flavour tetraquarks within the chiral constituent quark model 
of Glozman and Riska pT|-p^]. In this model, the hyperfme splitting in hadrons is due to 
the short-range part of the Goldstone boson exchange (GBE) interaction, instead of the one- 
gluon exchange (OGE) used in conventional models. In Ref. [p0|1 , we found that the GBE 
interaction strongly binds the heavy tetraquark system QQqq where Q = c or b ; this is at 
variance with the models based on OGE interaction where ccqq is unstable. Here we also find 
results somehow different from the previous literature. We show that the GBE interaction 
induces a short-range repulsion, of several hundreds of MeV in both P 5 and F\ systems. 
However, the GBE interaction also generates a long-range attraction due to its Yukawa-type 
potential. Together with a correlated two pseudoscalar-meson exchange, this will induce a 
long- and a medium-range attraction which could favour a deuteron-size (molecular type) 
bound state. In the following, we restrict ourselves to the study of compact objects. 

In the spirit of Jaffe's pioneering work here we raise the question whether or not the 
model of Glozman and Riska can accomodate a strongly bound, compact pentaquark. We 
shall compare our results with those based on OGE, as mentioned above. 

The model of Glozman and Riska reproduces well the light-quark baryon spectrum 
[2l|-p3| and some extension to heavy baryons has also been proposed ||24|| : in particular, 



it gives the correct ordering of positive and negative parity states in all parts of the consid- 
ered spectrum. The detailed form of the Hamiltonian used here will be given below. First 
it is useful to consider a simplified GBE interaction of the form (where radial dependence is 
neglected) : 

V^-C^Xf.Xf^.aj, (1) 

i<j 

where Af (F — 1, 2, 8) are the quark-flavour Gell-Mann matrices (with an implicit sum- 
mation over F) and a are the spin matrices. The minus sign of the interaction (|Tj) is related 
to the sign of the short-range part of the GBE interaction, crucial for the hyperfme splitting 



in baryon spectroscopy. This feature of the short-range part of the GBE interaction is dis- 



cussed at length by Glozman and Riska E1J . A typical order of magnitude for the constant 



C x is about 30 MeV. 



Here, we study pentaquarks with strangeness S — 0, —1, —2 and —3. As in Ref. PU| , we 
neglect the interaction due to meson exchange between a light quark and a heavy antiquark. 
Then the GBE interaction in a pentaquark q 4 Q reduces to the interaction between the light 
quarks q = u, d or s. 

The subsystem of four light quarks must be in a colour 3 state, or alternatively [211]c, 
in order to give rise to a colourless pentaquark. To calculate the GBE matrix elements, we 
assume that all light quarks are identical and that the ground-state orbital wave function is 
symmetric under permutation of light quarks. The corresponding Young diagram notation 
is [4]o- Then all combined spin [f] s and flavour [f] F symmetries, allowed by Pauli principle, 
form the states indicated in the first column of Table |. They have been obtained from inner 
product rules |25[ applied to the permutation group S , 4 . In column 2, the corresponding 
values of the spin S are indicated for each [f] s representation of SUs(2). The total angular 
momentum of the pentaquark is J = S + sq where sq = 1/2. The isospin I associated 
to a given SUp(3) representation is given in columns 3-6 for four distinct contents of the 
four-quark subsystem with strangeness S = 0,-1, —2 and —3. The last column indicates 
the expectation value of the operator (|I|) which can be conveniently calculated from the 
formula: 

(V x ) = N(N - 10) + + 1) + 2C F + AC C (2) 

o 

where N = A, Cc = 4/3 and Cp is the SUf(3) Casimir operator eigenvalue. One can see 
that the states in column 1 are given in an increasing order of (V^.). Thus |[22]5[211]_p) 
is most favourable (negative) for the structure uudsQ (I = 1/2, J = or 1) or udssQ 
(I = 0, J = or 1), |[31] S [22] F ) for uuddQ (I = 0, J = 1/2 or 3/2) and |[22] S [31] F ) for 
usssQ (I = 1/2, J = 1/2). The column (V x ) also suggests that other q 4 Q systems, with 1.1 
different from the above values, would be heavier, thus there is no reason to present results 



for such obviously unstable states. 



The GBE Hamiltonian has the form 22 



tf = £™*+£j^ - f^r + £ + S ^) , (3) 

with the linear confining interaction : 

V con{ (r ij ) = ~Xt-X c j Cr ij , (4) 
and the spin-spin component of the GBE interaction in its SUf(3) form : 

KF=1 

+ E ^(ry)Af Af + ^(ryJAjAj + fyfaW) *, • B h (5) 



F=4 



with A = y 2/3 1, where 1 is the 3x3 unit matrix. The interaction (H) contains j — it, K,r] 



or 77' meson-exchanges and the form of V^r%j) is given explicitly in Ref. [22] as the sum of 
two distinct contributions: a Yukawa-type potential containing the mass of the exchanged 
meson and a short-range contribution of opposite sign, the role of which is crucial in baryon 
spectroscopy. For a given meson 7, the meson exchange potential is : 

V"/( r ) = 7^" To {Q(r-r )fj, =a 3 exp(-a 2 (r - r ) 2 )}, 7 = 7r, K, 77, 77'. (6) 

47T YlmiVTij ' r ^Jtx 

For the Hamiltonian we use the parameters of Ref. P2| . These are : 

9l ± = 9% L = 9j<q = _ 6? flV, = L 

47T 47T 47T 47T 

r = 0.43 fm, a = 2.91 fm" 1 , C = 0.474 fm -2 , m ud = 340 MeV, 



/V = 139 MeV, n v = 547 MeV, /v = 958 MeV, fx K = 495 MeV. (7) 

They provide a very satisfactory description of low-lying nonstrange baryons, extended to 
strange baryons in Ref. |23| where a dynamical three-body calculation is performed as well. 
The latter reference gives m s = 0.440 GeV. 
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First, we calculate variationally the masses of the threshold hadrons using the Hamil- 
tonian (|3|)-(|7p. The meson mass is obtained from a trial wave function of type $ oc 
exp [— ar 2 /2] with r — f i — r 2 by varying the parameter a. The heavy-quark mass is 
adjusted to reproduce the experimental average mass M = (M + 3M*)/4 of M = D or B 
mesons. This gives m c = 1.35 GeV and = 4.66 GeV. With all quark masses now fixed, 
we estimate the baryon masses with a trial wave function of type $ oc exp [—a(p 2 + A 2 )/2] 
where p = r 2 — r*3 and A = (2r*i — r* 2 — r?) / v^- The results for baryons and mesons are exhib- 
ited in Table ||, together with the experimental values. The theoretical values lie at most 
~ 50 MeV above the experiment. The baryons masses, obtained from the same Hamiltonian 
but in a more precise Fadeev type calculations |23| are also indicated. 

For the pentaquarks, it is useful to consider the following system of internal Jacobi 
coordinates: 



(8) 



x — n - r 2 , 

y = (fx +f 2 - 2f 3 )/ v / 3, 

z = (fi + f 2 + f 3 - 3f 4 )/\/6, 

t = (ri + f 2 + r 3 + r 4 - 4f 5 )/v / 10. 

Then we assume a variational wave function of the form: 

* = f -V /4 ( -V ^xp h^(x 2 + + ^ 2 ) - ^t 2 ]. (9) 



7T/ \ 7T 



The expectation value of the kinetic energy T involves the average of the inverse light masses 
defined as: 

l/m a + 3/m b for q a ql 
= { (10) 
2/m a + 2/m b for g 2 g 2 , 



for the light subsystem. With 



^- = jU-!-, (ii) 

/i 2 4/xi m Q 



the internal kinetic energy of the q 4 Q system reads 



For the calculation of the confinement potential energy, we first need to evaluate the colour 
operator matrix elements : 

(0 C > = (Ef|> = Wc-4g, (13) 

i<j 

where C q = 4/3. For the colour state [211]c one has Cq = 4/3. Hence < 0° >[2ii]= —2. 
As there are 6 pairs in flT3| ) one gets 

<f.f> = -l/3, Z ' J = g ' (14) 

for each pair. For the antiquark the Casimir operator average is Cq = 4/3. Then for the 
light-heavy system, in a colour-singlet state with C c = 0, the colour average operator is 

(t y = {C ° C ~ Cc ~ C $ )/2 = _4/3 ' (15) 

i=l 

from which each quark shares 1/4. Hence 

<y-|> = "1/3, i = q,j = Q. (16) 

This proves that the colour matrix element is —1/3 for all pairs. With the trial wave 
function (Q) the space part of the confining interaction can be calculated exactly. Then the 
total confining energy is 

The matrix elements of the spin-flavour operators entering the GBE interaction (|5]) have 
been calculated by using the fractional parentage technique described in Ref . [|5|] . Accord- 
ingly, we first determine the explicit form of the spin and flavour states of a given symmetry 
[f]s and [f] F respectively. Then we rewrite each wave function as a sum of products of the 
wave function of the first two and of the last two particles. These are all either symmetric 
or antisymmetric two-body states, so that Eqs. (3.3) of Ref. [^TJ can be straightforwardly 



applied. In Table |T|, we exhibit (V x ) of Eq. (|5|) for the cases of interest. The upper index uu 
or us in the orbital two-body matrix element V designates the nature of the quark masses 
in the product miirtj in Eq. With the wave function (|9[), V can be explicitely calculated 
as: 



V* 3 



g 2 1 4 a 3 / 2 j> 2 



4vr 11 in, in J n 1 / 2 [2 a 



exp [-r (fi + ar )] - 



a: 



2 exp(a 2 rl) (a + a 2 ) 3 / 2 



4a 2 r 



a + a 2 ) -Jtx 



+ 2 + 



exp 



/i 2 /(4a) 



Erfc 



4a 4 r 2 \ 
a + a 2 J 

ji + 2 a r 
2v^ 



exp 



/x 3 a- 3 / 2 



4a 2 r 
a + a 2 

/4 



Erfc 



-a 2 r 



\J(a + a 2 ) 



(18) 



In Table pi , we present the energy E of each q^Q system under discussion, obtained by 
minimizing (H) with respect to a and b. Compared to the lowest threshold mass Et given 
in the last column of Table |TJ, one can see that the pentaquark energy lies several hundred 
MeV above the corresponding threshold, in all cases. One can also notice that the lowering 
of E — Et when replacing the quark c by b is small, of the order of 10 MeV. 

Among the systems studied here, the pentaquarks with S = —2 appear to be the most 
"favourable", i.e., have the lowest E — Et, although the value of E of a uudsQ system is 
lower than that of a udssQ system, due to a larger GBE attraction, as discussed above. The 
presence of two strange quarks, instead of one, has a negligible effect. Thus the lowering of 
E — Et for udssQ is essentially due to an increase in Et when passing from the iV + sQ 
threshold to the A + sQ threshold. 



In this first simple application of the chiral constituent quark model ||2T|— P4|| to the 
pentaquark spectrum, we did not find any deeply bound state. We however believe that 
the present results are not completely discouraging for the reasons which follow. The energy 
of strange pentaquarks is located above the threshold only roughly 1/2 from the amount 
found for hexaquarks with one heavy flavoured quark |2B| . The values found for E depend 
on the approximations involved. A lowering of E can be achieved by an increase of the 
basis vectors. A possibility would be to consider orbital states of type [22] containing two 
quanta of excitation, like in the nucleon-nucleon case [[^J. The lowest allowed value of V x , 



S 



Eq.(|), is -16 C x as for | [4] [22] S [211] F ) of Table |. Therefore an interference of the [4]o 
state considered here with the [22] o state is possible. 

Moreover the quark-antiquark interaction was entirely neglected here. Assuming that 
the instanton-induced interaction is very important for mesons (see e.g. [0), one could 
expect that an additional attractive contribution of this type would contribute towards 
stability of pentaquarks with respect to the meson-baryon threshold. As mentioned in the 
introduction, the GBE interaction through its Yukawa potential tail can also generate a 
medium- and long-range interaction between clusters. The structure of Eq. (J|) suggests 
that, in a frozen configuration, a strong repulsion at short distances implies a long-range 
attraction. Hence our conclusion on the non-existence of deeply-bound compact pentaquarks 
encourages pursuing theoretical investigation of deuteron-like baryon-meson systems, based 
on a dynamical approach of a few-body system, as for example the resonating group method. 
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TABLES 

TABLE I. The allowed |[/],s[/].f) states of colour 3 for the subsystem of pentaquarks with 
strangeness 0, —1, —2 and —3. The spin S and isospin / are indicated in each case |29|. The 
quantity (V x ) is the corresponding expectation value of the operator ([j]) in units of C x . 



MsWf) s I (V x ) 

uuddQ uudsQ udssQ usssQ 

|[22] 5 [211] F ) / 1/2 / -16 

|[31] s [211]f) 1 / 1/2 / -40/3 

|[31]s[2%) 10 1/21/ -28/3 

|[22] s [31]f) 1 1/2,3/2 0,1 1/2 -8 

|[4] 5 [3L» 2 1 1/2,3/2 0,1 1/2 

|[31]s[4] F > 1 2 3/2 1 1/2 8/3 
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TABLE II. Masses of hadrons required to calculate the threshold energy Et = baryon + 
meson. For D and B we indicate M = (M + 3M*)/4. For the baryons, a comparison is shown 
with the Fadeev calculations of Ref. |23|] 



Hadron 




Mass (GeV) 






variational 


Kef. J2J 


experiment 


AT 

IN 


U.y (\) 




0.939 


A 

A 


1 ICC 


1 1 OC 

1.13b 


i 1 1^ 
1.11b 


s 


1.235 


1.180 


1.192 




1.377 


1.348 


1.318 


D 


2.008 




1.973 


D s 


2.087 




2.076 


B 


5.302 




5.313 


B s 


5.379 




5.375 
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TABLE III. Energies of pentaquarks with J = 1/2 and 5 = 0, —1, —2 and —3. In column 3, 



(7 = 7r, 7], K or rf) designates the quark-quark matrix elements of the interaction 



Pentaquark I 



(V x ) 
(Eq.(D) 




Lowest threshold 



uuddc 



uuddb 



uudsc 



uudsb 



udssc 



udssb 



usssc 



usssb 



1/2 



1/2 



WV uu _ - 4/3 V3 



uu 
V' 



1/2 7 yuu _ 7 / 9 yuu _ U /Q V uu + 



22/3V% s + 22 /9V™ - 22/ 



28/3 V% s + 28/9 V™ - 28/9V3* 



V 2 22/3^ + 26/9^-26/9^ 



us _ 
V' ' 



20/9V„ ss - 10/9V3 3 
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